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Abstract 


The paper describes a set of results 
from the first phase of an RAE/NASA colla- 
borative programme on rotorcraft system 
identification that has the main objective 
of improving prediction methods. Plight 
measurements collected at RAE Bedford on 
an experimental Puma helicopter are 
reviewed and some notable characteristics 
highlighted. Following a brief review of 
previous work in rotorcraft system identi- 
fication, the results of state estimation 
and model structure estimation processes 
applied to the Puma data are presented. 

The results, which were obtained using 
NASA developed software, are compared with 
theoretical predictions of roll, yaw and 
pitching moment derivatives for a 6 degree 
of freedom model structure. Anomalies 
reported in other investigations have 
reappeared in this study. The theoretical 
methods used are described in the Appendix 
where a framework for reduced order model- 
ling is outlined. 


Notation 


A 

B 

CQi 


Crf 

^YFN 

F 


g 


XX 

’ yy’ 

L , 

L 

etc 


P 

M , 

M 

etc 


w 


%lc 

Nv» etc 

NrI 


P J qj r 


R 

r2 


-state matrix ( equation ( 7 ) ) 
-control matrix (equation ( 7 ) ) 
-rotor induced torque coef- 
ficient (equation (10)) 
-rotor thrust coefficient 
( equation ( 3) ) 

-fin sideforce coefficient 
(Fig 5 ) 

-force and moment vector 
( equation ( 7 ) ) 

-gravitational acceleration 
-height of aircraft eg below 
rotor (ft,m) (equation (7)) 
-moments of inertia of air- 
craft in roll, pitch and yaw 
(slugs ft^, kg m2) 

-rolling moment concise deri- 
vatives (normalised by Ixx^ 
-pitching moment derivatives 
(normalised by lyy) 

-pitching moment derivative 
wrt flapping (normalised by 


■-yy 
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-yawing moment derivatives 
(normalised by Izz) 

-yawing moment from main rotor 
(normalised bylzz) 

-roll, pitch and yaw rates 
(rad/s) 

-rotor radius (ft,m) 

-(multiple correlation coef- 
ficient ) 2 
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-solidity 

-time to half amplitude (s) 

-trim values of aircraft forward 
and normal velocity components 
(ft/s, m/s) 

-velocity perturbations along 
body x,y and z axes (ft/s, m/s) 
-sideways and vertical velocity 
components (ft/s, m/s) 

-control vector (equation ( 7 )) 
-state vector (equation (7)) 
-sideslip angle (rad) 

-rotor longitudinal and lateral 
disc tilts respectively (equa- 
tion (1) ) 

-Look number (equation ( 1 )) 
-relative damping 
-longitudinal and lateral stick 
position {% aft, std) 

-pedal position left) 
-longitudinal and lateral cyclic 
pitch angles (equation (1)) 
-aircraft pitch angle 
-system eigenvalue 
-rotor downwash (normalised) 
velocity (equation (10)) 

-flap frequency ratio 
-normal velocity component at 
rotor disc (equation (10)) 

-air density 
-relative air density 
-aircraft roll angle 
-rotorspeed 

-dutch roll frequency (equation 

( 9 )) 


Introduction 


A collaborative programme between 
the Royal Aircraft Establishment and the 
NASA - Ames Research Center to develop and 
exchange information on rotorcraft system 
identification is now underway with the 
main objective of improving prediction 
methods. In the first exercise of this 
kind at RAE, flight tests have been made 
on a Puma helicopter and the results have 
been analysed at NASA - Ames, using soft- 
ware developed in support of the Rotor 
Systems Research Aircraft programme. 

This paper reports on the first phase of 
this activity. The paper begins by des- 
cribing the aircraft and data processing 
system at RAE Bedford and attempts a pre- 
liminary interpretation of selected flight 
data based on comparison with theory. 
Aspects of rotorcraft system identifica- 
tion are then reviewed and the techniques 
currently in use are described. 
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Derivatives, estimated by these techniques 
for the Puma at a nominal 100 kn trim con- 
dition, are presented and compared with 
results predicted by a linear, 6 degree of 
freedom theoretical model developed at 
RAEI. Some of the major anomalies are dis- 
cussed. Finally the use of reduced order 
linear models for handling studies is 
reviewed, and approximations to the natural 
modes discussed in some detail. 


Flight Mechanics Investigations with 
PumaXW2tjl at RAE Bedford 

The Aircraft and Data Processing System 

Puma XW 241 (Pig 1) is a multi- 
purpose experimental aircraft operated and 
managed by the Helicopter Section of the 
Flight Research Division at RAE Bedford. 
Since its procurement in 1974, the Puma 
has been used in a variety of research pro- 
grammes and has recently been fitted with 
a digital Pulse Code Modulated (PCM) 
recording system; the results described in 
this paper were obtained with this equip- 
ment. A block diagram highlighting the 
features of the current data acquisition 
and processing system is shown in Fig 2. 

The instrumentation system includes 
three packs of linear accelerometers of 
both ac pickoff and force feedback type, 
two packs of rate and attitude gyros and a 
heading gyro. Air data sensors give air- 
speed measurement from a conventional 
pitot static and incidence and sideslip 
angles from vanes located on a nose boom. 
All four control displacements, the three 
swash plate jacks and the tail rotor pitch 
jack are sensed by potentiometers. In 
addition, one of the four blades is cur- 
rently instrumented to measure flap, lag 
and feather angles at the bearings. Samp- 
ling rates vary from 32 to 256 per second, 
the lower rate restricted to slowly vary- 
ing quantities, eg airspeed. The current 
overall rate for the multiplexed system is 
8024 words (12 bit) per second but the 
pattern of signals in each data field is 
flexible . 

During the Spring of 1981 a trial 
was conducted with the object of data 
gathering for a system identification 
exercise. Trims and responses to pilot 
shaped control inputs over a range of 
flight conditions were recorded. Data 
tapes from one of these flights, designa- 
ted 325, were sent to the Ames Research 
Centre for analysis using NASA software 
during the first round of the RAE/NASA 
collaborative programme in Rotorcraft 
Flight Mechanics. Results described in 
this paper are drawn from this exercise 
and used in subsequent interpretation, 
analysis and comparison with theory. 


Data pertaining to the trim condi- 
tions for flight 325 are given in Table 1. 
The nominal trim IAS was 100 kn from which 
a series of shaped control inputs were made 
and the ensuing response measured and re- 
corded. The inputs included steps, doub- 
lets and "3211" multisteps and generally 
produced repeatable response patterns; 
atmospheric conditions at the test points 
were very smooth . 

Time to recovery varied with input type and 
size but typically step inputs gave about 
10 s of data and 'return-to-trim' multi- 
steps gave about 20-30 s duration. 

Typical results from cyclic and pedal in- 
puts are reproduced in Figs 3 and 4. Before 
introducing the system identification pro- 
cess a preliminary interpretation will be 
attempted on these results together with a 
limited comparison with theoretical 
predictions . 

Preliminary Interpretation of Flight 
Results and Comparison with Theory 

An initial assessment of the flight 
results and comparison with theory provides 
a suitable background to judge the effect- 
iveness of both the theory itself and the 
system identification method, described 
later. The theoretical results presented 
in this paper are based on a 6 degree of 
freedom linearisation of the nonlinear 
simulation model described in Ref 1. 

During the development of this nonlinear 
model, a degree of validation was achieved 
through flight/theory comparisons of Puma 
data; these results are reported in Refs 
1 and 2 . 

The short term pitch and roll rate 
response to longitudinal and lateral cyclic 
steps inputs respectively are shown in 
Pig 3. The comparisons with theory are 
encouraging and indicate that the principal 
damping, control sensitivity and 'static 
stability' parameters ought to be predicted 
with reasonable accuracy by theory. There 
is some evidence that the 'initial' angular 
acceleration is sharper in theory but this 
is to be expected with a quasi-steady rotor 
model. This phenomenon and the difficul- 
ties it can present to a derivative estima- 
tion process are discussed in more detail 
in the next section. 

Turning now to the flight results 
shown in Pig 4, we see the coupled response 
to a doublet input applied to the pedal. 

This set of data reveals the presence of a 
lightly damped lateral/directional oscilla- 
tion that was very apparent to the pilot and 
observer performing these tests. This 
oscillatory 'dutch roll' mode is also pre- 
dicted by theory and a comparison of the 
measured and predicted characteristics is 
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shown in Table 2. It can be seen that the 
lateral/direotional characteristics are 
roughly in agreement whereas the longitu- 
dinal couplings, particularly the phase 
relationships 5 are significantly different. 
Having established possible basic ingred- 
ients of this mode it is useful to examine 
the theory to determine the origins. The 
derivative matrix used in the analysis is 
given in Table 3 and the system eigenval- 
ues in Table . Included in Table are 
the longitudinal and lateral subsystem 
eigenvalues derived with coupling terms 
set to zero. Conventional fixed-wing air- 
craft terminology is used to describe the 
different modes. The frequency of the 
dutch roll oscillation is determined pri- 
marily by the directional stability, -ie 


“dutch roll ~ 

An important effect displayed in 
Table 4 is the considerable reduction in 
damping of this mode when the coupling 
terms are included. The damping is rather 
low even for the lateral subsystem and 
this can largely be attributed to the low 
value of the yawing moment derivative Ny . 
This derivative contributes to the mode 
damping when the effective centre of the 
oscillation is a condition of non-zero 
sideslip. Wind tunnel data^ for the Puma 
fin in sideslip, as used in the theoreti- 
cal model, is shown in Fig 5. The side- 
force is strongly nonlinear with sideslip, 
but over the range -5° < B < +8°, practi- 
cally no lift is produced by the fin. It 
is believed that this effect is due to the 
suction on the rear of the 'lower surface' 
at small angles of attack; a characteris- 
tic of thick aerofoil sections with large 
trailing edge angles^. This aspect is 
discussed further in Ref 5, in relation to 
the Puma and how design improvements 
obviate the effect on the AS 332 Super 
Puma . 

Returning to the coupled system, a 
sensitivity analysis reveals that the main 
coupling derivatives affecting the dutch 
roll damping are Mp and . For a 

perfectly governed rotorspeed the deriva- 
tive reflects the torque changes pro- 

duced by the powerplant, following rotor- 
speed variations due to incidence perturba- 
tions. The effect of this derivative on 
the dutch roll and longitudinal short per- 
iod eigenvalues is illustrated in Fig 6. 

Nw is destabilising for the dutch roll 
mode and the roots are close to the longi- 
tudinal and lateral subset approximations 
when is zero. Pig 7 illustrates how 

the phase relationship between w and r 
(incidence and yaw rate) perturbations vary 
with for the dutch roll mode. They 

are seen to converge as increases and 

clearly when they are in phase, torque 


reductions impose a positive yawing moment 
on the fuselage when the yaw rate is a maxi- 
mum, giving a negative damping contribu- 
tion. This effect is discussed in more 
quantitative terms in the Appendix. The 
above description of the increased insta- 
bility of the lateral oscillation is de- 
pendent on the correct phase relationship 
between longitudinal and lateral variables. 
As highlighted earlier this is the area of 
most serious discrepancy between flight and 
theory. It should perhaps be emphasised 
that engine/rotor dynamics are not Included 
in the theoretical model but that there is 
evidence® that including this degree of 
freedom can reduce the effective yaw 
damping. Included in Fig A is the main 
rotor torque variation (via shaft strain 
gauge measurement) during the oscillation. 
Although fairly noisy, the fluctuations 
are seen to be phased relative to yaw rate, 
so that as the aircraft is yawing to star- 
board, the engine is applying a clockwise 
torque (minimum) to the fuselage. This is 
in the same sense as predicted by theory 
and it may be surmised that the phase shift 
relative to pitch rate and vane oscillations 
is accounted for in the engine dynamics. 

As a final note in this exploratory 
section some measurements were made on the 
Puma following a pedal doublet input with 
the pitch axis stability augmentation 
switched on, to give some indication of 
the effect of reducing the longitudinal 
coupling. A comparison of yaw and pitch 
vane responses with and without augmenta- 
tion is shown in Fig 8. Clearly the damp- 
ing of the oscillation has increased with 
the reduced longitudinal coupling. In 
fact, the time to half amplitude has now 
reduced to about one cycle, ie approxi- 
mately the same as predicted by the lateral 
subset calculation. The frequency is prac- 
tically unchanged. Longitudinal coupling 
obviously does play an important role in 
the damping of this oscillation. 

R otorcraft System Identification 

The state of the art in the field of 
system identification has largely been 
developed in the fixed wing community, the 
most recent comprehensive review of which 
can be found in Ref 7. There have been 
several attempts to apply the various 
techniques to helicopters, a review of 
which also appears in Ref 7. Before des- 
cribing the methods used in the RAE/NASA 
collaboration it is perhaps worth making a 
few observations on these past efforts and 
highlighting some of the lessons learned. 

Observations on Previous Work 

The ground rules for helicopter 
system identification were laid down in the 
pioneering work by Molusis in the early 
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seventies^ . In his papers, Molusis 
emphasised the need for long data records 
and combined manoeuvres in order to provide 
sufficient information for reliable esti- 
mation. He also stressed the importance 
of providing satisfactory initial esti- 
mates of derivatives in maximum likelihood 
algorithms by using an optimal filter/ 
smoother in conjunction with a least 
squares estimator. Of particular rele- 
vance to the helicopter problem were 
Molusis' s observations on the effects of 
rotor/fuselage coupling on derivative esti- 
mation. These were clearly demonstrated 
by attempting an identification of a 6 
degree of freedom (dof) model from simula- 
tion data that included rotor flapping 
modes!®. The identified derivatives were 
substantially different from quasi-steady 
predictions even though time histories 
showed reasonable correlation. Including 
the three principal multi-blade coordinates 
in the assumed model and performing the 
reduction to a 6 dof model after the iden- 
tification, resulted in very good agree- 
ment with quasi-steady values. The con- 
clusion was not that 6 dof quasi-steady 
theoretical models were necessarily inade- 
quate but that profound difficulties could 
be expected in trying to estimate these 
from flight measurements of fuselage 
motion alone. Molusis suggested an empiri- 
cal correction based on computed theoreti- 
cal data but his own results showed that 
this was not very satisfactory!®. A 
perusal of the simulation results in 
Ref 10 reveals that perhaps the most ser- 
ious discrepancy is the underestimation of 
the primary rate damping derivatives by 
the 6 dof identification. These usually 
produce dominant effects about all axes 
and relatively simple theory should be 
fairly reliable. A gross underestimation 
of this effect will obviously lead to a 
further corruption of secondary effects in 
each equation when the time histories are 
forced to agree. A possible cause of this 
anomaly can be demonstrated by a fairly 
simple example. 


Consider a hovering rotor with zero flap- 
ping hinge offset and centre of gravity 
(eg) on the rotor shaft. The first order 
flapping equations in multi blade coordi- 
nates e2.c ^Is > that represent 

longitudinal and lateral disc tilt respect- 
ively, can be written as 
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where y is the rotor Look number, SJ the 
rotor speed, p and q the body roll and 

pitch rates and e ,0 the longitudinal 
IS Xc 


and lateral cyclic pitch respectively. 
Non-uniform inflow effects have been neg- 
lected for this example. Constraining the 
body to rotate only in pitch and neglecting 
forward speed effects, the above equations 
can be augmented by the body pitching 
moment equation 
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lyy is the pitching moment of inertia and 
the pitching moment derivative Mg2c 

be written in terms of the rotor thrust 
coefficient C^ , 

Mgi^ = - hj^p(nR)^irR^C^ . (3) 


Here hp is the eg distance below the 
rotor, p the air density and R the 
rotor radius . For present purposes we 
have assumed in equation (3) that the 
rotor thrust remains normal to the disc 
during pitching motions. The quasi-steady 
form of equation (2) can be written in the 
form 

q = M^q + 01 s (4) 
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From equation (1) the quasi-steady 
derivatives can be written as 
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To demonstrate the phenomenon in question we 
now attempt to estimate the derivatives Mq 
and ■’'50]_g in a model structure given by 

equation (4) from data generated by equations 
(1) and (2), following a step input in e^g , 

Configuration data used for the numerical 
study are y=8.0, a = 27.8 rad/s and 
= - 6 . 36 . Using a least squares 

estimator, the effect of data length on the 
estimated derivatives is shown in Fig 9. A 
sampling interval of 0.04 s was used for 
this case. The effect is dramatic, particu- 
larly on the estimated damping Mq . Using 
only a short data length a positive Mq is 
predicted and even after 2 s the estimation 
is -Still only 70% of the quasi-steady value. 
For a 2 s data run, the effect of sampling 
interval on the estimation is shown in 
Fig 10. Increasing the sampling interval 
to 0.2 s results in an accurate estimation 
of the quasi-steady derivatives. Obviously 
the model structure given by equation (4) is 
inadequate for portraying the short term 
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dynamic effects of the flapping equation 
(1), which in this case produces a mode 
with frequency 5-5 rad/s and time to half 
amplitude of about 0.1 s. A comparison of 
pitch rate response using the full equa- 
tions and the quasi-steady representation 
is illustrated in Fig 11. Clearly the 
angular acceleration in the former case is 
zero at the origin and builds up to the 
quasi-steady value after about 0.2 s. 

Using data in this interval for estimating 
quasi-steady derivatives tends, therefore, 
to underestimate the damping in the motion 
to the extent that for, very short data 
lengths, an unstable system is predicted 
as shown in Fig 9 . 

When analysing flight data one would, 
of course, use a much longer data record 
than used in the previous example and the 
step input is rather limited in its excita- 
tion spectrum. However, a multi-step input 
would typically extend over about 30% of 
the data record and the contamination 
effects described above might persist well 
into the record. 

Several other attempts to identify 
6 dof helicopter models have resulted in 
underestimation of primary damping deriva- 
tives relative to the quasi-steady values 
(Refs 11-13) . All these references report 
using advanced statistical methods for 
estimating derivatives but the phenomenon 
is still present in the results. The 
obvious solution to the problems high- 
lighted by Molusls is to measure individual 
blade motions and to use higher order model 
structures in the estimation process. There 
will obviously be flight cases when such a 
solution is mandatory to the Identification 
of certain phenomena, particularly when 6 
dof theoretical models are known or sus- 
pected to be inadequate. For a wide range 
of conditions within the normal flight 
envelope, however, it is hoped that 6 dof 
models can be used to predict handling 
characteristics adequately and it is worth- 
while exploring other potential solutions 
to the data processing problems before 
resorting to the measurement of blade 
motions. One approach is to use an input 
signal with a frequency content that 
excludes the rotor modes. A possible 
shortcoming of this method is that the data 
could be starved of information necessary 
to predict some of the derivatives, par- 
ticularly for hingeless rotor helicopters 
where the lower frequency rotor modes are 
not far removed from some of the body 
modes. Nevertheless, the use of smooth 
inputs rather than sharp steps or ramps 
has an obvious appeal in this context. A 
similar and perhaps complementary solution 
involves filtering the data to exclude the 
'higher' frequency content but again this 
could mean coming down to frequencies as 


low as 2 Hz. Some exploratory work was 
done on this theme by one of the present 
authors in Ref 14, where the estimation 
process was transferred to the frequency 
domain. The data can then, be filtered by 
excluding data beyond the frequency range 
of interest. This idea will now be exploi- 
ted for the analysis of the Puma data. 

Current Methods 

The techniques used in this paper are 
based on the system identification approach 
described in Ref 15, where the data analy- 
sis is divided into three main stages; 
state estimation, model structure estima- 
tion and parameter identification. The 
state estimation stage involves the recon- 
struction of all states and controls from 
the available sensor measurements using a 
Kalman filter/smoother^° . This process 
also eliminates the biases and scale factor 
errors from the measured data and reduces 
the level of measurement noise. The second 
stage utilises a stepwise least squares 
regression techniquel7 to derive an ade- 
quate model structure from a range of pos- 
sible candidates . The method uses certain 
statistical criteria to determine which 
states should be included in the model and 
estimates the corresponding coefficients. 
The third stage adopts these estimates as 
starting values for a maximum likelihood 
estimation process that generates unbiased, 
efficient, parameter estimates. 

The analysis of results described in 
the next section will be confined to the 
first two stages described above. The 
equations used in the state estimation pro- 
cess are essentially kinematic relation- 
ships. For the present data set the recon- 
structed states include the three velocity 
components and the translational and rota- 
tional accelerations of the aircraft centre 
of gravity. The external forces and moments 
are then directly related to the accelera- 
tion estimates and for the model structure 
estimation we choose the 6 dof formulation 
of the form 

F = Ax + Bu (7) 

P , X and u are the external force and 
moment vectors, vehicle state vector and 
control vector respectively. The regres- 
sion process currently used treats each 
equation separately, but rather than find 
the coefficients of this limited model that 
are valid for the entire frequency range of 
the data, it is more appropriate to deter- 
mine the best set for the frequency range 
of interest. This is readily accomplished 
by transforming the data to the frequency 
domain using a Past Fourier Transformation 
(PPT) and then truncating the data beyond 
the bandwidth of interest. If it is 
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assumed that the coefficients do not vary 
with frequency then the regression can be 
performed in the frequency domain with the 
vectors in equation (7) interpreted as the 
corresponding Fourier transforms. An 
additional advantage of estimation in the 
frequency domain is that the resolution 
can be improved by padding the data with 
zeros. This effectively forces the spec- 
tral lines to be more closely spaced with- 
out affecting the frequency content of the 
data. The higher resolution should result 
in more accurate estimation of the low 
frequency modes, which are difficult to 
identify from a limited duration data 
recordl^ . 


Model Structure Estimation with Puma Data 

The results presented in this sec- 
tion were produced at the Ames Research 
Center and forwarded to RAE Bedford for 
interpretation. The inherent communica- 
tion difficulties in this process limited 
the scope of this first phase of the col- 
laboration. The analysis software is cur- 
rently being implemented on a VAX computer 
at Bedford when these difficulties should 
be largely overcome. 

From the assortment of data collec- 
ted during Puma flight 325, three man- 
oeuvres were selected for analysis, namely 
doublet inputs in lateral cyclic, pedal 
and longitudinal cyclic. The data used 
consisted of measurements from one of the 
packs of Inertial instruments (accelero- 
meters, rate gyros, attitude gyros), the 
airspeed sensor and the pitch and yaw 
vanes. The instrument pack referred to 
included a 2-pole Butterworth filter 
(ojn ~ 10.6 Hz , ? = 0.73) that effectively 
removed the dominant 4/rev noise in the 
measurements. The results of the model 
structure estimation process reveal a 
marked variation in estimated coefficients, 
particularly the contributions from coup- 
ling effects. The reasons for this are 
not yet understood and we therefore choose 
to concentrate on a synthesis of the prim- 
ary moments for each input, 7e rolling 
moment for lateral cyclic doublet, yawing 
moment for pedal doublet and pitching 
moment for the longitudinal cyclic doub- 
let. The sample rate for all channels 
was 64/s, a rate that limited the record 
lengths to' about 15 s. This limitation 
will be removed in future analysis. The 
results from each manoeuvre will be dis- 
cussed in turn. 

Lateral Cyclic Doublet/Rolling Moment 
Synthesis 

Results from the filter/smoother 
process are shown in Fig 12, where it can 


be seen that only two cycles of the lat- 
eral oscillation are available. Two fur- 
ther points are worth noting; a bias has 
been detected and corrected for in the 
pitch rate gyro and the initial slope of 
the normal velocity component (w) indicates 
an initial climbing and decelerating flight 
condition. The effect of this unsteady 
initial condition has not been explored. 

In discussing the results of the model 
structure estimation we will refer to the 
multiple correlation coefficient (R) for 
the fit. In essence, the closer r 2 is to 
unity, the better the overall fit of the 
data in the frequency domain. Table 5 
summarises the estimated rolling moment 
coefficients (derivatives) for the cases 
studied, along with the current quasi- 
steady theoretical predictions. Deriva- 
tive estimates with a structure containing 
only lateral variables are compared for 
three frequency ranges, 0 to 0.5 Hz, 0 to 1 
Hz, 0 to 4 Hz. In all cases the high r 2 
values indicate a reasonably good fit with 
this limited model. The roll damping Lp 
is seen to increase by about 50% as the 
bandwidth is reduced but is still markedly 
lower than the theoretical prediction. The 
Fourier transform of the rolling moment is 
compared with the estimated fit in Pig 13, 
for the two lower frequency ranges. Most 
of the data is contained in the range 0 to 
0.5 Hz, and the dominant peak at 0.25 Hz 
corresponds to the lightly damped lateral 
oscillation. The fairly close fit at and 
above this frequency appears to be achieved 
at the expense of the fit at lower frequen- 
cies. Referring again to Table 5, it can 
he seen that R2 rises to 0.96 when the 
longitudinal variables w and q are 
added, with some associated modification to 
the lateral variable coefficients. Compar- 
ing the estimated and theoretically predic- 
ted derivatives in Table 5 we can see that 
the major anomalies are for the rate der- 
ivatives Lp and Lr ■ The magnitude of 
these effects seems to have been reversed 
for the flight results but it is difficult 
to rationalise the physical significance 
of this; indeed the result seems rather 
dubious. Further exploration into the 
details of the regression analysis is 
clearly required. 

P edal Doublet/Yawing Moment Synthesis 

The smoothed state estimates for this 
manoeuvre are illustrated in Fig 14 and 
the yawing moment derivatives derived from 
the lower frequency bandwidth of data are 
compared with theoretical predictions in 
Table 6. Results from two different model 
structures are included, one with lateral 
variables only and the other with the addi- 
tion of the longitudinal variables w and 
q . The multiple correlation coefficient 
is hardly affected by these additions. 
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It can be seen that the directional stiff- 
ness (Ny) and damping (N^) compare very 
well with theory but the control power 
(N^ ) is somewhat larger according to 

theory. The most striking differences are 
in the coupled rate derivatives Np and 


N, 


The small value of Np predicted by 


theory is the result of the cancellation 
of two larger effects from the main (nega- 
tive Np) and tail-rotor (positive Np) , 

It is interesting to note that the value 
of r 2 was raised to O. 9 I by the first two 
variables drawn into the regression, v 
and rip j the remaining variables then 
being added to account essentially for the 
out of phase component of the response. 

For the lightly damped yawing mode the sum 
of these contributions should obviously be 
small. Once again, a detailed breakdovrn 
of the in-phase and quadrature components 
of these 'damping' terms is required to 
aid further interpretation of the results. 


L q ngi tudinal Cyclic Doublet/Pitching 
Moment Synthesis 

The state estimates derived for this 
manoeuvre are shown in Pig 15, where the 
bias error in the pitch rate gyro is again 
apparent; also the normal and longitudinal 
velocity estimates again indicate an initial 
decelerating climb. The least squares 
estimates in the frequency domain are com- 
pared, in Table 7, with the quasi-steady 
theoretical predictions for the derivatives. 
Results are presented for the longitudinal 
model structure for the three frequency 
ranges. Both the pitch damping (Mq) and 
control power (M^iqg) are seen to increase 

as the bandwidth is reduced and both rise 
to nearly 80% of the theoretical predic- 
tions. The remarkable agreement for the 
speed stability derivative (M^) is somewhat 
overshadowed by the incidence stability 
derivative (M^j) comparison. The variable 
w actually just managed to become part of 
the fit, accounting for only the last few 
per cent in the value of R2 . It is pos- 
sible that the estimation of an unstable 

is related to the unsteady initial con- 
ditions for this manoeuvre, as the effects 
of the initial deceleration persist for the 
duration of the record as shown in Fig 15. 

It is intended to explore this topic 
further . 


The results from the three manoeuvres 
discussed above are both encouraging and 
perplexing. Unfortunately, time was not 
available to explore further the uncovered 
anomalies, for this paper. The underesti- 
mation of roll damping, relative to theory, 
a feature common to several earlier 
attempts at rotorcraft parameter estima- 
tion, is perhaps of greatest concern. The 
reduction effect of the lower frequency 
flapping modes, described in an earlier 
section, should be minimized by truncation 


in the frequency domain. The least squares 
regression analysis does of course have its 
shortcomings, particularly when noise is 
present in the dataTo^ ljut also when there 
is strong correlation between states. How- 
ever the technique is appealing in that it 
offers a simple and systematic approach to 
model structure estimation, and investiga- 
tions of the type described above will be 
pursued in the continuing collaborative 
programme. When fully exploited and under- 
stood, the results of the model structure 
estimation stage will then be used to 
initiate the more complex maximum likeli- 
hood process. 

The model structure estimation stage 
serves another useful purpose in the vali- 
dation of reduced order approximate models 
of helicopter flight mechanics. Simplified 
models, which still indicate trends 
accurately, have obvious advantages and 
various schemes are suggested in the 
Appendix for the 6 degree of freedom heli- 
copter as a framework for analytic model 
reduction. The results described are based 
on the theoretical predictions of Puma 
characteristics discussed earlier in the 
paper. Clearly, however, we have failed to 
validate the theoretical arguments put for- 
ward to explain the low damping of the 
lateral oscillation. In particular, the 
yawing moment derivatives N^, and Nq in 
Table 6 bear little resemblance to their 
theoretical counterparts, and, unfortunately 
for this study no estimate was obtained for 
the derivative Mp (another important 
effect in the theory) . There is no reason 
to believe that any major physical effect 
has been neglected in the "theoretical quasi- 
steady derivatives, but it is possible that 
the omission of the rotorspeed degree of 
freedom has distorted the model structure 
estimation. With this possibility in mind, 
future estimations will include this addi- 
tional state. 

Concluding Remarks 

Results from the first phase of an 
RAE/NASA collaborative programme in rotor- 
craft system identification have been des- 
cribed. Flight data from three manoeuvres 
with an RAE Puma helicopter have been pro- 
cessed by the state estimation and model 
structure estimation processes of a NASA 
system identification software package. 

A comparison of moment derivative estimates 
with theoretical predictions has been used 
as a guide to the likely accuracy involved. 
No firm conclusions can be drawn from this 
first exploration but a number of features 
are worth highlighting. 

(1) Encouraging results have been obtained 
by performing the regression analysis in 
the frequency domain. Reducing the band- 
width of data used resulted in an increase 
in roll and pitch damping and control power, 
as expected from time domain considerations. 
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butj in some cases, the results still fall 
well short of expected quasi-steady values. 

(2) Primary yawing moment derivatives 
estimated from a pedal doublet manoeuvre 
show good agreement with theory but the 
adverse yaw (Np) and coupling from longi- 
tudinal motion (Nv,,,Nq) show significant 
differences. The contribution of these 
coupling terms to the 'dutch roll' mode 
damping, so elegantly expressed by simple 
theory, cannot, therefore, be 
substantiated . 

(3) Roll damping (Lp) and control power 
(Lnic) estimated from a lateral cyclic 

doublet manoeuvre are considerably lower 
than theory predicts - a result in keeping 
with previous reports. This anomaly in the 
rolling moment is accompanied by a very 
high estimate for the derivative Lp and 
the temptation at this stage is to question 
the structure estimation, rather than 
theory, for these effects. 

(4) The closest agreement with quasi- 

steady predictions were obtained for pitch- 
ing moment derivatives from a longitudinal 
cyclic doublet manoeuvre, except for the 
static stability derivative which is 

estimated to be very small and positive by 
the least squares regression. 

(5) A theoretical framework for exploring 
reduced order model structures has been 
outlined in an Appendix to the paper. 

A coupled longitudinal/lateral fourth order 
system is required to describe the current 
theoretical predictions of the Puma 'dutch 
roll' oscillation. 

The collaborative programme is still 
at an early stage and more detailed 
investigations using the model structure 
estimation process are planned for differ- 
ent manoeuvres of longer duration and poss- 
ibly improved model structures. 

Appendix 

The Use of Reduced Order Theoretical Models 

The complex nature of helicopter 
flight mechanics make it a prime candidate 
for treatment as a sum of interacting sub- 
systems, with the attraction that phenomena 
may be described by considering a series of 
lower order problems and their interactions. 
Conditions under which this type of approx- 
imation is valid are often based on intui- 
tive reasoning but they can be formulated 
more precisely using notions of subsystem 
dynamic separation (eg widely separated 
characteristic times) and Interaction 
strength. For aircraft flight dynamics 
such opportunities often arise for des- 
cribing rigid body /aeroelastic interaction 
or in the description of the individual 
modes using familiar arrangements of motion 
variables, eg longitudinal short period mode 
made up of incidence and pitch rate excur- 
sions. Once again, the attractions of 


performing such a 'partitioning' analysis 
can often stem from the enhanced physical 
appreciation it inspires . In this 
Appendix we concentrate on the theoretical 
results discussed earlier and outline a 
framework for reduced order analysis. One 
can clearly imagine the 6 dof model as 
already being a reduced order approximation 
of a higher order system containing rotor 
as well as engine/fuel control system 
degrees of freedom. 

The method of analysis is described 
in Ref 19 where the concept of weak coup- 
ling is introduced and conditions of 
application are quantified. The technique 
has been applied to strongly controlled 
aircraft motions^O and more recently to 
describe the range of application of the 
longitudinal short period approximation for 

helicopters 21 . jn the present paper, a 
form of approximation for lateral/direc- 
tional motion is sought. For the Puma 
results described the strong coupling from 
longitudinal motion renders the search in 
vain. In the following, aspects of open 
loop stability characteristics only will be 
addressed . 

The reduction process is based on a 
partitioning of the system matrix into 
lower order subsystems that are weakly 
coupled. Details of the method can be 
found in any of the references cited above 
and will not be elaborated upon here. In 
order to achieve the correct partitioning 
for rigid body modes of motion we need, to 
introduce the vertical velocity wq and 
sideways velocity vq as new variables, 
replacing pitch attitude and yaw rate in 
the equations of motion. The system 
matrices, in partitioned form, for decoupled 
longitudinal and lateral motions are shown 
in Table 8, along with the approximating 
characteristic polynomials. One would not 
expect these formulae to give very accurate 
results in general but they do serve as a 
yardstick against which the effects of coup- 
ling terms can be measured. For the Puma 
derivative data given in Table 3 the weakly 
coupled approximate results are shown in 
Table 9. 

A comparison with the decoupled longi- 
tudinal/lateral results given in Table 4 
indicates again that the dutch roll damping 
is badly overpredicted. This poor compari- 
son will result when the oscillatory side- 
slip motion has a significant component of 
sideways motion superimposed on the side- 
slip due to yawing motion. The overpredic- 
tion of the phugoid damping can be attri- 
buted to a similar effect in the longi- 
tudinal plane. For the other three modes 
the approximations are clearly satisfactory. 
The importance of the directional stiffness 
Ny on the dutch roll damping is not, of 
course, predicted at all by this type of 
approximation. The root loci in Fig 16 
illustrates this effect; the two curves 
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shown are for the coupled longitudinal/ 
lateral system and for the lateral system 
alone. Both loci are approaching the 
asymptotic value given in Table 9 as Ny 
increases, ie dutch roll damping “(Nr +Yy). 
The point being made is that as directional 
stability increases the formulae given in 
Table 8 are not only becoming better 
approximations to the lateral subset damp- 
ing but that both are improving relative 
to the fully coupled result. In these con- 
ditions the other modes are also approxi- 
mated fairly accurately, eg spiral mode 
marginal stability is well predicted, 
although these results are not shown here. 

The results discussed above indicate 
that the degree of coupling between longi- 
tudinal and lateral motion can be strongly 
influenced by the directional stiffness or, 
in other words, by the frequency of the 
lateral oscillation. This result is some- 
what intuitive and for the present coupled 
system the coincidental similarity between 
the longitudinal short period and dutch 
roll frequencies is bound to result in 
strong longitudinal/lateral interactions. 

As discussed in an earlier section these 
interactions are brought about mainly 
through the derivatives Mp and ; the 

main rotor contribution to these deriva- 
tives being negative and positive res- 
pectively for 'clockwise' rotating rotors. 
The strong coupling present leads to a 
further reduction in dutch roll damping 
through the mechanism discussed earlier. 
This effect can be quantified by consider- 
ing a reduced order system made up of both 
the longitudinal and lateral 'fast' 
oscillatory modes. 

Assuming on the one hand that trans- 
lational velocity excursions (u,vo,Wo) are 
much slower and weakly coupled, and, on 
the other, that the roll subsidence 
approximation Interacts in a quasi-steady 
manner, the 4th order approximate system 
takes the form. 


w 


z 

Z + U 

0 



w 

q e 


q 


M 

M 

-M L /L 



w 

q 

p V p 

V 


0 

0 

0 

V 


-N U 

0 

UN + (L /L ) ( 

— 


L w e 


^ e V V p 


Here we have neglected all coupling terms 
except Mp and , but with the strong 

coupling remaining, no further reduction 
would be reliable. Examination of the 
Routhian for this 4th order system indi- 
cates that the oscillatory stability bound- 
ary is crossed when N^„ ~ 0.02 , which is 
in accordance with the results presented 
earlier in Pig 6. The eigenvalues for the 


above system are, 

^short period 
^dutch roll 


= - 1.0 ± 1.0791 

= - 0.114 ± l,115i 


giving the right order of damping reduction 
for the dutch roll mode. If we assume a 
neutrally stable oscillation exists for the 
above system then a good approximation to 
this damping decrement can be obtained by 
deriving the steady state frequency res- 
ponse of equation (8). For this case w 
is approximately l 80 ° out of phase with 
V (hence in phase with yaw rate) and the 
effective damping becomes. 


da"Ping(dutch roll) 

- («r ^ ^ ^ - NpU^) 


N M O^L 

w p e V 

^p“0 ^^w ’^q^ 


(9) 

“ 0.234 , 

hence 

R (X) - - 0.117 

e 

(where Uq is the dutch roll frequency) 

which agrees with the result given by 
equation (8) above. 


The derivative , as stated 

earlier, represents the quasi-steady torque 
variation produced by the engine, in res- 
ponse to rotorspeed variations . The 
validity of this implicit weak coupling 
assumption can only be assessed when an 
engine/rotorspeed control system model 
structure, representative of Puma, is 
incorporated. The derivative itself 

is produced mainly by the 'induced torque', 
written in coefficient form. 




- /q) 


( 10 ) 


M L /U L 
P r e p 


X 

g-NU),N + Y +(L/L)(g-NU) 
^ p e/ ’ r V r p p ej 


( 8 ) 


Crji is the rotor thrust coefficient and 
(fz ~ ^0^ upwash, normal to the disc. 

The variation in the semi-normalised yawing 
moment (yaw acceleration), from this source, 
with normal velocity w is shown in Fig 17. 
The variation is seen to be moderately non- 
linear particularly in the normal helicopter 
working state (w < 0). The derivative 
is seen to increase as autorotation is 
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approached “ ^0 ~ 0) but the effect 

will, of course, disappear once the engine 
disengages. It is known that the direc- 
tional stiffness can also reduce In this 
region due to fin shielding effects; the 
loss of 'dutch roll' damping for these 
flight conditions should, therefore, be 
fairly severe if the current predictions 
are correct. 
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Notoinal* lAS 
Relative density^ 
Rotorspeed, n 
Rotor torque 
A(7tf 




I 


XX 

A* 


yy 


I 


*A 

ZZ 


a 


100 kn 

0.818-0.837 

26.6-27.2 rad/s 

1 1500 ft lb (564 hp, 421 kW) 

13018-12169 lb (5904-S519 kg) 

6380 slug (8650 kg m^) 

25483 slug (34550 kg 

20283 slug (27500 kg m^) 


Cg location below rotor hub 
Cg location forward of rotor hub 
Rotor Lock nuober 

Solidityg s 


7.05 ft (2. 157 m) 
0.086 ft (0.026 m) 

7.86 

1.0516 

0.0917 


Mode 

Coupled system 

Decoupled long/lat 
subsystems 

Roll subsidence 

-2.242 

-2.209 

'Dutch roll' 

-0.0451 il.047i 

-0.193 ±l.079i 

Spiral 

-0.1 166 

-0.1194 

•Short period* 

-0.9054 ±I.I86i 

-0.7645 i0.9354i 

'Phugoid* 

-0.00833 ±0.17641 

-0.0168 ±0.2038i 


Table 4. Comparison of coupled and uncoupled 
longitudinal/lateral system eigenvalues — 
Puma, 1 00 kn 


* Actual IAS varied from 95-100 kn giving an EAS 105-110 kn. 
** Manufacturer's estimates (Ref 3). 


Table 1 . Puma flight 325 — nominal trim conditions 


Parameter 

Flight 

(approx) 

Theory 

Period 

(s) 

4.5 

6.0 


(a) 

14.0 

15.36 

|p/v| 


0.02 

0.019 

tp/v 


135° 

158° 

[t/v| 


0.012 

0.0096 

4r/v 


-95° 

-85° 

|q/v| 


0.006 

0.0027 

tqlv 


98° 

-48° 

|w/v| 


0.336 

0.354 

^w/v 


30° 

-107° 


Table 2. Comparison of measured and predicted 
'dutch roll' oscillation characteristics 


u V 5 » P f 

X -0.2650E-OI -0.93iOE-03 0.397QE0I O.I890E-02 O.IOIOEOI O.OOOOE 00 
2 -0.3190E-01 -0.7050E 00 0.1680E03 0.I450E-01 0.2100E 01 O.OOOOE 00 
X 0.2460E-02 -0.5570E-02 -0.83I0E 00 -O.SOOOE-03 -0.2070E 00 O.OOOOE 00 
T -0.8120E-02 0.2930E-02 0.1040E01 -O.I360EOO -^.3980EOI -O.I690EO3 
L -0.5520E-O2 -0.2730E-03 0.3390E-00 -0.2200E-0I -0.2050E 01 0.2940E 00 
N 0.2880E-02 0.I270E-0I -0.3280EOO 0.6050E-02 -fl,9290E-03 -0.5280E00 


Table 3. Theoretical estimates of Puma derivatives 



Flight estimates 
lateral variables only 

+ 

w and q 

Theory 


0-0.5 Hz 

0-1 Hz 

0-4 Hz 

0-0.5 Hz 



0.92 

0.91 

0,9 

0,96 


L 

V 

-0.013 

-0.013 

-0.013 

-0.0123 

-0.022 

L 

p 

-0.4 

-0.32 

-0.28 

-0.35 

-2.05 

L 

r 

1 .7 

1 .82 

1 .9 

3.03 

0,294 

L 

0.022 

0.021 

0.02 

0.019 

0.044 







L 

w 




0.0054 

-0.00027 

L 

q 




1 .39 

0.839 


Table 5. Rolling moment derivatives ~ estimates from 
flight and theory 



Flight 

(0-0 

Lateral 

variables 

only 

es timates 
.5 Hz) 

+ 

w and q 

Theory 


0.97 

0.98 


N 

V 

0.0053 

0.0071 

0.00605 

N 

P 

-0.6 

-0.53 

-0.0009 

N 

r 

-0.636 

-0.572 

-0.528 

N 

-0.027 

-0.028 

-0.043 

N 

V 


-0.003 

0.0127 

N 

q 


0.322 

-0.328 
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Table 6. Yawing moment derivatives — estimates from 
flight and theory 



Flight estimates 


Theory 



0-0.5 Hz 

0- 1 Hz 0-^ 

Hz 

r2 

0.802 

0.76 0.73 

M 

0.00243 

0.00243 0.00245 0.00246 

u 




M 

0.00092 

0.00085 0.00074 -0.0056 

w 




M 

-0.683 

-0.667 -0.648 -0.831 

q 




^nis 

0.0316 

0.0302 0.029 0.038 

Table 7. Pitching moment derivatives 

— estimates from 


flight and theory 



(Note mTz~* <■ U ) = M (Z + U ) - 

M Z ) 


W q 

e w q e 

q w 

Longitudinal subset 

X = (u,w ,w,q) ; w 

= „ - u e 



- 0 o 



X g/U 

u e 

X - g/U X - w ‘ 

w “ e q e 



Z 0 

Z Z 



u 

w q 



Z,, 0 

Z Z + U 




v; q e 



M 0 

M M 



u 

u q J 



2 Cphagoid) ; 


I (x - t?-M H (Z h U ) + (X - W )zli )\ 

1 V ^e/ \ ® ^ c u w / 

I “ U ) I 


* u ) 


K V - w ^ T _1 


M CZ + U ) 
w q u 


^ (short period); X - (Z^ * ^ 


(y = “ 0) 


subset 

X = {v 
- o 

v,Vsp); V 

o 

= V + U 

e 

0 

0 

Y 

V 


0 

0 

1 

0 

-N 

r 

-U N 
e V 

N + Y 
r V 

g - N U 
p e 

L /U 
r G 

L 

V 

-L /U 
r e 

1 

a. 


^ 7==5 g \ 

G (l N + 7^ L I 

\ p V vj 

,7 '' ' ("x " '‘r * r (f - "p)] 

• (■• • ^ (t ■ "I "• ■ 


X_ (roll subsidence) ; X « L 
B p 


Table 8. Approximate formulae for longitudinal and 
lateral eigenvalues 


Roll subsidence 
'Dutch roll' 
Spiral 

'Short period' 
'Phugoid' 


-2.05 

-0.346 ±1.1 16i 
-0. 1124 

-0.768 ±0.968i 
-0.0349 ±0.1969i 


Table 9. Approximate eigenvalues for longitudinal and 
lateral subsystems 








1^- 


Fig 1 Puma XW241 



Fig 2 Puma flight mechanics data acquisition and 
processing at RAE Bedford 
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Fig 3 Flight/theory corhparison of rate response to 
cyclic inputs 
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Fig S Puma fin sideforca coefficient vs sideslip angle 



Fig 6 Root loci for fast oscillatory modes with varying 
yawing moment due to incidence (Ny^) 
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Fig 4 Puma flight 325, Run 12 - response to pedal doublet 
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Fig 7 Variation of the phase angle in incidence and yaw 
rate excursions during 'dutch roll' oscillation with 
Nyy — theoretical predictions 
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Fig 8 Puma - yaw and pitch vane response to pedal doublet 



Fig 9 Variation of derivative estimates with record length 



Fig 10 Variation of derivative estimates with sampling interval 



Fig 1 1 Comparison of pitch rate response with and without 
flapping dynamics 





2.0 4 0 6.0 8.0 10. 0 Us) 12.0 

Fig 12 Comparison of measured (— ) and estimated ( — ) 
states (Kalman filter/smoother) for lateral cyclic 
doublet manoeuvre 
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Fig 13 Fourier transform of measured (— ) and estimated ( — ) 
rolling moment — lateral cyclic input 
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Fig 16 Root loci for 'dutch roll' modo with varying fjg •\-j Variation of 'induced' rotor torque reaction 

Ny — comparison of coupled and lateral/ with normal velocity perturbations 

directional subset modes 





